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Using general scaling arguments combined with mean-field theory we investigate the critical (T ~ 
T c ) and off-critical (T / T c ) behavior of the Casimir forces in fluid films of thickness L governed by 
dispersion forces and exposed to long-ranged substrate potentials which are taken to be equal on 
both sides of the film. We study the resulting effective force acting on the confining substrates as a 
function of T and of the chemical potential fi. We find that the total force is attractive both below 
and above T c . If, however, the direct substrate-substrate contribution is subtracted, the force is 
repulsive everywhere except near the bulk critical point (T c ,/x c ), where critical density fluctuations 
arise, or except at low temperatures and (L/a)(/3Afi) = O(l), with Afi = fj, — fj, c < and a the 
characteristic distance between the molecules of the fluid, i.e., in the capillary condensation regime. 
While near the critical point the maximal amplitude of the attractive force if of order of L~ d in the 
capillary condensation regime the force is much stronger with maximal amplitude decaying as L~ x . 
In the latter regime we observe that the long-ranged tails of the fluid-fluid and the substrate- 
fluid interactions further increase that amplitude in comparison with systems with short-range 
interactions only. Although in the critical region the system under consideration asymptotically 
belongs to the Ising universality class with short-ranged forces, we find deviations from the standard 
finite-size scaling for fln(f/£^) L even for i,L > where £(t = (T — T c )/T c -> ±0, = 
0) = £ol*l _I/ ! i s the bulk correlation length. In this regime the dominant finite-size contributions 
to the free energy and to the force stem from the long-ranged algebraically decaying tails of the 
interactions; they are not exponentially small in L, as it is the case there in systems governed by 
purely short-ranged interactions, but exhibit a power law decay in L. Essential deviations from the 
standard finite-size scaling behavior are observed also within the finite-size critical region L/£ = O(l) 
for films with thicknesses L < L C rit, where -Lcrit = Co = (16|s|)"' / ' 3 , with v and /3 as the standard bulk 
critical exponents and with s = O(l) as the dimensionless parameter that characterizes the relative 
strength of the long-ranged tail of the substrate-fluid over the fluid-fluid interaction. We present the 
modified finite-size scaling pertinent for such a case and analyze in detail the finite-size behavior in 
this region. The standard finite-size scaling behavior is recovered only for L ^> L CT it- 
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INTRODUCTION 



A. General properties 



Since its first prediction in 1948 by Hendrik Casimir 
[H, the effect named after him has raised significant the- 
oretical and experimental interest and has been found 
to occur in numerous manifestations (see, e.g., Refs. 
@, IE II IE IE 0, B ©])■ Originally, Casimir consid- 
ered vacuum fluctuations of the electromagnetic field be- 
tween two parallel metal plates which restrict and mod- 
ify the fluctuation spectrum leading to a dependence 
of the energy on the distance L between the plates. 
This so-called classical (actually quantum mechanical) 
Casimir effect, which for decades had been considered 
as a theoretical curiosity, in recent years has turned 
into a subject of intensive research, not at least trig- 
gered by possible applications in micromechanical de- 



vices, and has been already verified with an impressive 
experimental accuracy [TcL fill . E2] • Another manifesta- 
tion of the Casimir effect, which has a significant impact 
towards a different scientific direction (see , e - g., Refs. 

@, 0, EE EI EE EE EE EE EE EE HI HI EE m HE EE) , 

is the so-called statistical-mechanical (thermodynamic) 
Casimir effect. In a fluid close to a continuous phase 
transition at a critical point T c large fluctuations of the 
order parameter occur. If, as before, the fluid is confined 
by parallel plates at a distance L, and is in contact with a 
particle reservoir with a chemical potential ji, the grand 
canonical potential fl ex (T, /i, L) of the fluid in excess to 
its bulk value ALw^i^T, fi) depends on L so that one 
can define the effective force / between the plates per 
cross sectional area A and per fcgTas 



f(T,n,L) = -0 



du>ex(T, n, L) 
dL 



(1.1) 
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where (3 = l/(k B T), v ex (T, (x,L) = uj(T,(j,,L) - 
LuJbui^T, fi) = Sl ex (T, fi, L)/A is the excess grand canon- 
ical potential per cross sectional area A, fi(T, fi, L) — 
Au>(T, /i, L) is the total grand canonical potential, and 
Ubuik(T, n) is the density of the bulk grand canonical 
potential. Besides temperature, chemical potential, and 
film thickness the force also depends on which boundary 
conditions the surfaces impose on the system. The or- 
der near the surfaces can be either reduced or - which is 
the generic case for liquids confined by solid substrates - 
increased due to effective surface fields generated by the 
confinement. The latter case is known as (+, +) bound- 
ary conditions (for a more precise definition see below). 
For this case the schematic phase diagram of a fluid film 
with thickness L is shown in Fig. [TJ 
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FIG. 1: The schematic phase diagram of a d-dimensional slab 
with (+, +) boundary conditions for various thicknesses L and 
surfaces which prefer the liquid phase. The bent gas-liquid 
coexistence curves correspond to capillary condensation tran- 
sitions for L — Li and L = L2 with Li < L2 (compare, c.f., 
Fig. 11(b)). Away from the critical region the shift in the 
phase boundary relative to the bulk coexistence line A/i = 
is proportional to L~ , while within the critical region it is 
proportional to L~ A ^ where A and v are the standard bulk 
critical exponents. The lines of first-order phase transition 
end at (d — l)-dimensional critical points T c (Li) with coor- 
dinates {T c ,Li, A/ic^J, i = 1,2, the positions of which vary 
with L and depend on the presence and on the strengths of 
the fluid-fluid and the substrate-fluid interactions. For large 
L these points are located close to the bulk critical point 
T c with coordinates (T c ,A[i = 0): T C} l — T c ~ L~ x l" and 
A/i Cj l ~ L~ A ^" . Since the fluctuations in systems of reduced 
size are stronger, one typically has T Ci Li < T c ; the fact that 
A^ic^i < expresses the preferences of the identical walls. 

One should keep in mind that the force /(T, /i, L) is a 
definition dependent quantity in the sense that it depends 
on how one defines the thickness of the film. A natural 
choice would be to take the distance between the planes 
defined by the positions of the nuclei of the top layer of 
each substrate. However, there are certainly other possi- 
ble definitions, which will differ by a microscopic length. 
This implies that a quantitative comparison between ex- 
perimental data and theory is only possible if the data 



are accompanied by a precise definition of what L is. Up 
to now seemingly there is no awareness of this issue yet. 

Close to the critical point, the critical Casimir inter- 
actions are proportional to IzbTc and therefore the inter- 
action between the plates can become rather strong in a 
system with high T c such as, e.g., in classical binary liquid 
mixtures. However, in such systems the direct dispersion 
forces between the fluid particles and between the fluid 
particles and the substrate particles play also an impor- 
tant role. The contributions of the dispersion forces to 
the total effective force can be distinguished from that of 
the critical Casimir forces by their temperature depen- 
dence, because the leading temperature dependence of 
the former does not exhibit a singularity. For this rea- 
son both in theoretical analyses and in interpretations of 
experimental data the contributions due to the disper- 
sion forces are usually treated separately and, also for 
the critical region, are simply added as a regular back- 
ground contribution to the total force (see, e.g., Refs. 

[3, [13, HI). 

B. Scaling with dispersion forces 

The present investigation aims at studying in detail 
the actual interplay between the dispersion and the crit- 
ical (fluctuation induced) forces in simple fluids bounded 
by strongly adsorbing identical walls. In describing fluid 
systems near criticality one often resorts to the Ising 
model as a representative of the corresponding universal- 
ity class, for which only nearest neighbor interactions are 
considered. Instead, here we consider long-ranged pair 
interactions between the fluid particles decaying asymp- 
totically ~ J 1 r~ 6 for distances r between each other 
and long-ranged substrate potentials ~ J l,s z~ 3 acting 
on fluid particles at a distance z from the flat surface of 
a semi-infinite substrate. (Here we do not embark into 
the extended microscopic description of confined binary 
fluid mixtures, which are expected to belong to the same 
universality class as the one-component fluids considered 
here explicitly.) Such systems also belong to the Ising 
universality class characterized by short-ranged forces 
[27I ] . This implies that the corresponding critical expo- 
nents and leading finite-size dependences are expected to 
be the same as for systems with short-ranged forces. The 
long-ranged part of the interactions leads to a regular 
contribution to the force, as considered and expected be- 
fore, and in addition, as shown recently [13, HI, HI, ES] , to 
a serious modification of the finite-size behavior of its sin- 
gular contribution. Indeed, in the case of (i-dimensional 
systems, in which the interaction decays asymptotically 
with the distance ~ r - d -<* (with a > 2), the effective 
total force /(T, /i, L) between the plates can be cast (see 
Eq. CD])) into the form 0, Hi, Ell, [11 

f(T,»,L)~ (1.2) 
L- d X\LH u Lfa, (L/aoV"' I, (L/toV^ s, (L/fo) - " g u ] 
+(3{a-l)H A {T,ii)L-°§- d . 
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FIG. 2: (a) Temperature dependence at A^t = of the total Casimir force normalized by its value at the bulk critical point in a 
system with L = 500 layers for different strengths A and s of the long-ranged part of the fluid-fluid and fluid-substrate potentials, 
respectively; L and £t are measured in units of the lattice constant of the lattice model. The total force depends significantly 
on the strengths of both potentials. Nonetheless it turns out that the critical part of the force exhibits universality and scaling 
behavior. For small L, however, important nonuniversal contributions due to the long-ranged parts of the aforementioned 
interactions emerge which cannot be neglected. Independent of the values of s and A the force extremum always occurs at 
positive t — (T — T c )/T c , i.e., above the bulk critical temperature. This is due to the stabilizing effect of both the boundary 
conditions and of the long-ranged tails of the interactions on the order of the system so that the strongest fluctuations are 
expected to occur slightly above T c and on the vapor side of the bulk coexistence curve. Since the Casimir force is generated 
by these fluctuations one expects and indeed observes that the force is strongest for such thermodynamic parameters, (b) Same 
as in (a) for the field dependence A/i at T = T c . The force extremum always occurs at negative values of A/i, i.e., on the gas 
side of the bulk coexistence curve. 



Here X is dimensionless, universal scaling function, £t is 
the bulk correlation length £(t -> ±0, A^ = 0) = 
at bulk coexistence fi = [i c and for t = (T—T c )/T c — > ±0, 
while Z„(t = 0,A^ -> 0) = Co, M l/3c Afi\-"/ A is the bulk 
correlation length at the critical temperature T = T c 
with j3 c = (fee^c) . For T > T c and A// = one has 
Co = ; while for T < T c and Afi = one has £o = £o~ 
with the ratio £o~/£o~ being universal. £o,^ is the same 
for Afi — > +0 and A/j, — ► —0. The second term in Eq. 
(|1.2|) stems from the free energy contribution HaL*" -1 ' , 
wheTeH A (T,n)=A t (T,iJ,)+Ai(T,fj,)+Ai !S (T) with so- 
called Hamaker constants At, Ai, and Ai, s (see below for 
more details). In addition, to is the standard Wegner's 
correction-to-scaling exponent for short-ranged systems, 
while uji = a — (2 — r\) and uo s = a — (d + 2 — rj)/2 
are the correction to scaling exponents due to the long- 
ranged tails of the fluid-fluid and substrate-fluid interac- 
tions, respectively. Further L-dependent contributions to 
the force such as next-to-leading order contributions to 
the Hamaker terms or higher order corrections to scaling 
are neglected because they are smaller than those cap- 
tured in Eq. Q1.2)) . Here we consider a d-dimensional 
system characterized by an interaction given as the sum 
of a short-ranged and a long-ranged component with the 
latter decaying according to a power law r~ d ~ cr , with 
2 < a < 4. We intend to vary the ratio A of the strengths 
of the long-ranged and the short-ranged contributions. 
By varying A we can quantitatively probe the importance 
of the long-ranged tails. One might envisage potential 
experiments in colloidal systems which allow for a ded- 
icated tailoring of the form of the effective interactions 
between colloidal particles. We recall that, independent 



of A, for 2 < a < 4 the system still belongs to the corre- 
sponding short-ranged universality class, i.e., the critical 
exponents do not depend on a. Such long-ranged in- 
teractions are called subleading long-ranged interactions 
0, H§| • In Eq. {15} Tj is the standard critical exponent 
characterizing the decay of the bulk two-point correlation 
function at the critical temperature, is the (dimen- 
sionless) scaling field associated with the Wegner-type 
corrections, while I and s are dimensionless nonuniversal 
coupling constants: I is proportional to the strength J 1 of 
the long-ranged tail in the fluid-fluid interaction whereas 
s is proportional to the contrast between the substrate 
potential and the fluid-fluid interaction integrated over 
a half-space (see below). For the "genuine" non-retarded 
van der Waals interaction, which governs nonpolar fluids, 
one has d = a = 3 and thus uj s > 0, lui > 0. We will fur- 
ther suppose that the positivity of uj s and wi is fulfilled 
for all values of d and a considered in the following. Ac- 
cordingly, for L large enough, one can expand the scaling 
function X in Eq. lfL2]l : 

/ ~ L- d {x SI [L/£ u L/^] 

+ (L a XtlL/tuLfo] 
+ (L/Co)- u "l Af [L/£ t ,L/^] 
+ guXZlLf&L/tp]} 
+P(a-l)H A (T,^L-^- d . (1.3) 

If Eq. I|1.3p is valid the scaling function X SI , which 
originates from the short-ranged interactions and de- 
scribes the well studied (short-ranged) critical behav- 
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ior, provides the leading behavior of the force near the 
bulk critical point (£/£< = 0, — 0). Here the no- 
tion "near the critical point" means L/£ t = 0(1) and 
L/£ M = 0(1), which defines the "critical region of the fi- 
nite system". There X^J, Xj r , and AT' 1 ' represent only cor- 
rections to the leading L dependence. In particular, the 
value X sr (0) = (d— 1) A Qi b is related to the Casimir am- 
plitude A a b, which depends only on the bulk universality 
class and the boundary conditions a, b at the two confin- 
ing surfaces, i.e., the surface universality classes. There 
is considerable knowledge about the Casimir amplitudes, 
such as their exact values in d = 2 [33L H3, Ha] and in 
d = 3 for the spherical model [l5| , e-expansion results for 
d = 3 [TUmi], and mean field values [la]. In the following 
we shall consider only the case of (+, +) boundary condi- 
tions corresponding to strong equal surface fields acting 
at both surfaces, i.e., the case of strong adsorption of 
the fluid on the boundaries of the bounding substrate. 
While the scaling function X s1 ' decays exponentially for 
L/£ ^S> 1, in this regime the other scaling functions Xf 
and X] r in Eq. I|1.3p decay according to a power law. For 
this reason these functions, which formally stem from 
corrections to scaling due to the subleading long-ranged 
tail in the fluid-fluid interaction and in the substrate po- 
tential, respectively, lead to leading finite-size contribu- 
tions in the singular behavior of the force in the regime 
L/£ 3> 1. In the case of periodic boundary condition - in 
which only the functions Xj r and X^ are present - the 
above statement has been verified for the spherical model 
(which, for periodic boundary conditions, represents the 
limit n — > oo of the 0(n) models). We note that although 
the scaling function X£ r characterizes the short-ranged 
universality class, the scaling field g u incorporates, in 
general, also contributions due to the long-ranged tails 
of the interaction. We refer the interested reader to Ref. 
[H| where explicit results about the above mentioned 
mixing of the corrections due to the long-ranged forces 
and the Wegner-type corrections have been reported. In 
the following we are only interested in the leading In- 
dependence of the force. Therefore we shall not discuss 
in further detail those contributions which are due to the 
corrections to scaling ruled b y th e Wegner exponent u> 
[36| ; they produce corrections f37j, [38| both near the crit- 
ical point, where they are of the order of , as well as 
away, where they decay exponentially. 



C. Relevance-irrelevance criterion 

We briefly comment on the conditions which justify an 
expansion of the type presented in Eq. f| 1 .3[) . The re- 
quirements 2 — T) — a < and (d + 2 — 77) /2 — a < 
are obvious and, as already mentioned, normally they 
are satisfied in any realistic system for which d = a = 3 
and 77 <C 1 (e.g., for the 3d Ising model 77 ~ 0.034 36]). 
One important additional condition arises (see Ref. 3Sjj), 
however, from the fact that we consider not semi-infinite 
systems but systems which are finite in one dimension 



and thus exposed to power law long-ranged substrate- 
fluid potentials from both sides. Accordingly these po- 
tentials act everywhere in the finite system; at the center 
of the film their sum is minimal but not zero. The effect 
of this value of the total external potential in the center 
can be interpreted as if the system has a nonzero effec- 
tive bulk excess chemical potential A/i G ff despite the fact 
that the actual bulk system might be at bulk coexistence 
curve A/a = 0. Taking into account the contributions 
from both surfaces in terms of the notations already in- 
troduced one has A/x e fl = 2s [L/(2£ )P <T - Since the ex- 
cess chemical potential scales as A/iL A ^ u one finds that 
in a film the finite-size contributions due to the long- 
ranged substrate-fluid potentials will be negligible in the 
critical region if 



2| S | [V(2£o)P [L/^f /u « 1, 



i.e., 



2 <h 



A/u—a 



< 1. 



(1.4) 



(1.5) 



The sign of s is chosen such that s > (< 0) corre- 
sponds to attractive (repulsive) walls, i.e., walls prefer- 
ring the liquid (gas) phase of the fluid. A more detailed 
discussion of that point will be given below where we 
identify s within the framework of a mean-field model. 
Due to standard relations between critical exponents 
one has Ajv = (d + 2 — r/)/2, so that relation (|1.5p is 
consistent with Eq. l|1.3p . On the other hand due to 
Ajv — a = d — cr — (3/v and with d — a for realistic 
systems the condition 1|1.5|1 turns into 



*H-ii 



-(3/v 



< 1. 



(1.6) 



With £ typically of the order of 3 A, (3 ~ 0.329, v ~ 
0.631 (3d Ising model), and for a = 3 one finds 

L » L crit = Co (2° +1 \ S \Y /0 * 612 H 1 - 918 A. (1.7) 



Equation l|1.7p represents a relevance-irrelevance crite- 
rion for van der Waals type substrate potentials in a 
film of thickness L: if L ^> L cr ;t the leading behavior 
of the force within the critical region will be that for a 
system with short-ranged interactions, while for thinner 
films L < L cr it the effects due to the van der Waals sub- 
strate potentials are expected to be relevant. Since s is 
typically not very small (see below), L cr a turns out to 
be surprisingly large. As discussed later for a variety of 
systems |s| £ [1,2]. However, for some systems such as 
3 He or 4 He films near their bulk liquid-gas critical point 
and confined by Au substrate s can be as large as 4 [39| . 
Thus for most fluid systems, for which measurements 
of the Casimir force has been performed, L was much 
smaller than L C rit- We stress, however, that until now 
no Casimir force measurements for systems with (+, +) 
boundary conditions near liquid-gas critical points have 
been reported. Therefore we are unable to confront our 
predictions explicitly with available experimental data. 
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Finally, we note that if for a particular choice of the 
fluid and of the surrounding substrate s happens to be 
small, also £ C rit, as given by Eq. (|1.7p . will be small. 
Keeping in mind that in Ising-like three-dimensional sys- 
tems (with d = a = 3 and rj g= : 0.034) the correction-to- 
scaling exponents u> ~ 0.81 [3a], uji ~ 1.03 and oj s ~ 0.52 
are numerically not very different from each other, the de- 
termination of L cr it should in principle in this case take 
into account also the values of the scaling fields g u and 
I in addition to that one of s, i.e., Eq. (|1.7p should be 
modified. However, in such a case an expansion of the 
type given in Eq. (|1.3p will be valid for any reasonably 
large L with the contributions due to the Wegner type 
corrections as well as due to the long-ranged part of the 
interaction competing with each other and representing 
corrections to the leading behavior of the force which is 
given by the scaling function X ST . 



D. Total Casimir force 

Figure [2] presents the typical behavior of the total (in- 
cluding the contributions from the regular part of the free 
energy but neglecting the direct substrate-substrate in- 
teraction) Casimir force as a function of the temperature 
or the excess chemical potential, respectively. The force 
is normalized by its value at the bulk critical point. The 
precise meaning of the parameters s and A, as well as the 
definition of the model within which the curves have been 
calculated will be given below. Larger values of A and s 
correspond to stronger long-ranged tails of the fluid-fluid 
and substrate-fluid interactions, respectively. The case 
(A = 0, s = 0) corresponds to truly short-ranged inter- 
actions. We expect these curves to resemble potential 
experimental data for the force occurring in non-polar 
fluid films bounded by substrates which strongly prefer 
the liquid phase of the fluid. (Experimentally the contri- 
bution of the direct substrate-substrate interaction can 
be separated off by a control experiment with the fluid 
replaced by vacuum; therefore in the following we do not 
consider this direct contribution.) Despite of the spread 
of the curves for different values of s and A it will turn 
out that a description in terms of scaling functions is pos- 
sible if the data are suitably analyzed, provided that for 
a given value of s the film thickness L is sufficiently large 
(see Eq. (T7|). 

It is our aim to determine the leading L-dependence of 
the Casimir force as a function of the temperature, the 
chemical potential, and of the parameters s and A, and to 
clarify the relevance of these parameters for its behavior 
in a fluid film bounded by strongly adsorbing, identical 
walls. To this end and in view of the large parameter 
space, starting from a continuum density functional we 
describe the system in terms of a lattice gas model which 
is expected to capture the essential features of actual 
fluids as far as the critical behavior is concerned. For 
L ^> L crit the terms entering into Eq. I|1.3|) are the most 
important ones (36l. l37l. l38| . because for the systems under 



discussion here they provide the dominant singular finite- 
size behavior of the free energy and of other thermody- 
namic variables. Whenever s^0, which is the generic or 
so-called normal case for actual fluids confined by solid 
substrates, the functions X sr and X l J will determine the 
leading critical finite-size behavior of the Casimir force. 
[We note that if the fluid film of interest is a wetting 
fluid on a solid substrate or an interfacial wetting film 
in binary liquid mixtures, one or two confining specta- 
tor phases are fluids, too. Apart from subtleties such as 
the occurrence of capillary waves at confining fluid inter- 
faces, our present approach is expected to cover the im- 
portant aspects also of these more complicated systems. 
The different preferences of the walls, however, lead to a 
situation which is not described by (+, +) boundary con- 
ditions and requires, therefore, a separate consideration.] 
Due to the lack of the symmetry s — » — s in the presence 
of (+, +) boundary conditions, a term in the free energy, 
which is - to lowest order - linear in s, is indeed ex- 
pected. So-called neutral walls, considered theoretically 
(see, e.g., Ref. [40] and references cited therein), corre- 
spond to s = 0. In this special case X sr and Xj T will 
determine the dominant finite-size behavior of the fluid 
system. For such a case X sr and X] T must be calculated 
for Dirichlet-Dirichlet boundary conditions. Recently a 
similar system of a fluid near a weakly attractive wall has 
been considered by Monte Carlo methods [41] ; but in this 
study the interaction potential between the fluid particles 
has been truncated thereby rendering effectively a short- 
ranged potential. Finally, we emphasize again that for 
film thickness L < L CTl t an expansion of the type given 
in Eq. f| 1 . 3fl is not possible so that for these systems the 
van der Waals interaction will be important throughout 
the critical region, including the point (T = T C , Afi = 0), 
i.e., in this case even the determination of the Casimir 
amplitudes is affected. 



E. Outline 



The article is arranged such that in Sec.[ll]we present 
the lattice gas model with van der Waals type long- 
ranged interactions. In Sec.[ni] we present our results 
for the free energy and the force by focusing in Sub- 
sec. |lIIA]on our results for the force near criticality, while 
in Subsec. BlIBI we discuss the behavior of the system off 
criticality, including the capillary condensation regime. 
Finally in Sec. lIVI we summarize our findings and com- 
ment on their experimental relevance. As far as possible 
we compare our results with previous findings for fluids 
governed by or exposed to dispersion forces. Important 
technical details are presented in Appendices [Bj and 

ej 
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II. THE MODEL 

Within the density functional approach for inhomoge- 
neous fluids, which in practice is mean-field-like in char- 
acter, the grand canonical functional 0[p(r)] of a fluid 
has to be minimized with respect to the local number 
density p(r) [13, S3| (in view of our focus on critical phe- 
nomena we refrain from using more sophisticated versions 
of density functionals) : 



n[p(r)} 



/hsWf)) 



(2.1) 



1 



p(r)w(r - r')p(r')d 3 r d 3 r' 



+ J V{z)p{v)d 3 r- p J p{r)d 3 r. 

The fluid is confined between two parallel flat plates at a 
distance L which exert a substrate potential V(z) with z 
as the normal distance from one wall. For an individual 
wall Viz — > oo ) ~ z~ a with a = 3 for a genuine van der 
Waals interaction; p is the chemical potential and fns(p) 
is the bulk free energy density of a hard-sphere system 
acting as a reference system; w(r) is the fluid potential 
which is suitably regularized at a molecular distance, i.e., 
it exhibits a negative, finite, plateau around |r| = and 
tends to zero according to a power law for large values 
of r = |r|. In Eq. l|2.ip the integrals run over the slab 
volume. 

In the spirit of focusing on the essential features of 
confined critical phenomena, one can further simplify the 
continuum functional in Eq. (|2.1j) by replacing it by its 
lattice version. This resembles the approach taken, e.g., 
by Fisher and Nakanishi |45| in their mean-field study 
of the same system but governed by short-ranged forces. 
The grand potential functional for this lattice gas system 



fi[p(r)] 



k B Tj2{p(r)^[p(r)} 
rec 

+ [l-p(r)]ln[l-p(r)]} 
+ \ £ p(r)w(r - r')p(r') 



■J2[V(z)-p]p(r), 



(2.2) 



rec 



where C is a simple cubic lattice in the region < z < L 
occupied by the fluid. Here and in the following all length 
scales are taken in units of the lattice constant a of the 
order of a molecular diameter (and thus are dimension- 
less) so that the particle density p(r) is dimensionless and 
varies within the range [0,1]. In Eq. l|2.2p the terms in 
curly brackets correspond to the entropic contributions, 
while in an obvious way the other terms are directly re- 
lated to the interactions present in the system. 

The variation of Eq. (|2.2p with respect to p(r) leads 



to the equation of state for the equilibrium density p* (r) 
V(r)-1 (2.3) 



tanh 



2 2^ w{ - r 



r')p*(r') + £(p-V(z)) 



The advantage of this equation is that it lends itself to 
solving it numerically by iterative procedures. For a 
given geometry and surface potential V(z) the solution 
determines the equilibrium order-parameter profile p* (r) 
in the system. Inserting this profile into Eq. l|2.2p renders 
the grand canonical potential of the system. In order to 
avoid the double sum in Eq. (|2.2p , which is inconvenient 
for the numerical treatment, from Eq. (|2.3p one can eas- 
ily derive the relation 



k B Tj2p*(r) arctanh[2p*(r) - 1] (2.4) 
rec 

lj2(»-V(zW(r)-± ]T w(r-v')p*(r)p*(r'), 



rec 



r.r'eC 



which yields for Eq. (|2.2 

n[p*(r)] = £ 

rec 

+ [1 - p*(r)] In [1 - p*(r)] - p*(r) arctanh[2p*(r) - 1]} 



k B T{p*(r)ln[p*(v)} 



±[p-V(z)]p*(r) 



(2.5) 



Note that here p*(r) is no longer a free functional vari- 
able, with respect to which one has to minimize Eq. (|2.5p . 
but the solution of Eq. (|23| . 

Denoting <j>*(r) = 2p*(r) — 1 and Ap = p — p c , where 
p c — \ J2 r ' w ( r ~ r ')> t ne equation of state (|2.3p can be 
rewritten in the standard form 



4>* (r) = tanh 



/?£ J ( r > r ') 



AV(z)) 



(2-6) 

where J(r — r') = — w(r — r')/4. The bulk properties of 
the model are well known (see, e.g., Refs. [40, [43] and ref- 
erences therein) . We recall that the order parameter 0* of 
the system has a critical value (j>* — which corresponds 
to p c — 1/2 so that 4>* = 2(p* — p c ). The bulk critical 
point of the model is given by (/3 = (3 C = [J2 r ^( r )]~\ P — 
p c = — 2J2r J( r )) with the sum running over the whole 
lattice. Within the mean-field approximation the critical 
exponents for the order parameter and the compressibil- 
ity are = 1/2 and 7 = 1, respectively. The effective 
surface potential AV(z) in Eq. (|2.6p is given by (see Ap- 
pendix [Aj in particular Eq. (j A8|l as an analogue) 

AV(z) = 5v s [(z + l)- a + (L+1 — z)- a ] , (2.7) 
1 < z < L — 1, where contributions of the order of z~ IJ ~ 1 , 



etc. have been neglected, 



8v s 



' 1 + 

, 2 , 



-( Ps J Ls - pcf) 



(2.? 
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is a (T- and /^-independent) constant, 

J(r) = Jl {5{\v\) + 5{\v\ l)} + J l /{l + \v\ d+ °) e{\v\-l), 

(2.9) 

is a proper lattice version of —w(r)/A as the interaction 
energy between the fluid particles, and 

J Ls (r) = J l s f 5(\r\ - 1) + J l >*/\r\ d +° 9(\r\ - 1) (2.10) 

is the one between a fluid particle and a substrate par- 
ticle (here S(x) is the discrete delta function [48| and 
6(x) is the Heaviside step function with the convention 
9(0) = 0); p s is the number density of the substrate par- 
ticles in units of a~ d . Note that Eq. (|2.9[) incorporates 
the lattice version of the regularization of J(r) at r = 
as implied by the original off-lattice density functional 
approach. Here we take into account on equal footing 
the two contributions to the effective interface potential 
stemming from the missing fluid particles and from 
the substrate particles. From Eq. I|2.6|l - (|2.10p one can 
identify the dimensionless coupling constant 



s = -^[3 5v s 



(2.11) 



(see Eq. (j 1 .3jl : s > corresponds to walls preferring the 
liquid phase of the fluid). The case s = 0, which will be 
considered later, corresponds to p s J ' s = p c J l . Note that 
the effective potential Sv s is formed by the difference of 
the relative strength of the substrate-fluid interaction for 
a substrate with density p s and that one of the fluid-fluid 
interaction for a fluid with a density p c . In Eq. I|2.7|l the 
restriction z > 1 holds because we consider the layers 
closest to the substrate to be completely occupied by the 
liquid phase of the fluid (which implies that we consider 
the strong adsorption limit), i.e., p(0) — p(L) — 1, which 
is achieved by taking the limit Jj,f — > oo; therefore the 
actual values of AV(0) = AV(L) will play no role. In 
order to preserve the monotonic behavior of w(r) as a 
function of the distance r between the particles, in Eq. 
(EH) we have to require that Jj r > J 1 /{I + 2^ d+a ^ 2 ). 
In terms of 6 the functional (12.21) turns into 



fi[#r)] = k B Tj2\ 
, 1 - <Hr) 



1 + 4>{r) 



In 



1 + 0(r) 



In 



1 - <p(r) 



-i^[A/z-AF(*Mr) (2.12) 

re-C 



r,r'£C 



where 



2 ^ 



Ap-AV(z)-J2 JM 



r'eC 



(2.13) 



does not depend on 4> and therefore is a regular back- 
ground term which carries a L-dependence and thus 
shows up in the force. An expression similar to the one 
in Eq. (|2.5p . which avoids the double sum and thus is 
more convenient for numerical procedures, can be also 
obtained. With the identifications </>(r) «-» m(r) and 
|[A/i-Ay(z)] «-» h(z) one can rewrite the above expres- 
sion for f2[(/>(r)] as a functional Afi[m(r)] = (Q — 0, xeg ), 
which describes a magnetic system at temperature T and 
in the presence of an external local and spatially varying 
magnetic field h(z). In the remainder we shall use this 
mutual correspondence between the fluid and magnetic 
systems in order to make contact with existing theoretical 
results for any of them. We shall call the force calculated 
from using fi[^(r)] the total Casimir force, while that part 
of the force calculated with the regular background term 
subtracted will be termed the critical Casimir force. This 
procedure corresponds to the analysis of experimental 
results if one subtracts from the value of the force mea- 
sured around T c the asymptote obtained by measuring 
this force well above T c . 

In accordance with Eq. (|1.3p , for the finite-size behav- 
ior of the excess grand canonical potential per unit area 
A of a liquid film in the case when both confining surfaces 
prefer strongly the liquid phase, one expects 

We»(t, Mi L\d, a) ~ + < 2 (2.14) 
+k B TL- (d - 1) x 

x a [L/z u (pa^)(l/^) a/i/ , (m )-"% (m )-"*s] 

+ [Ai(T, p) + A ltS (T, p) + A S (T)} L-^-V, 

where er"^ and cr" s 2 are the non-singular parts of the sur- 
face tensions at the surfaces 1 and 2, respectively (the 
singular parts are incorporated in Xn). Here 



MT,P) 



47r (rf~D/2 r( i±g) 



J l p 2 b (T 7 p), (2.15) 



with pb being the bulk fluid density at the given T and 
p, represents that part of the Hamaker constant which 
is generated by the long-ranged part of the fluid-fluid 
interactions, i.e., the dispersion interaction [49, 50], A s < 
is the part due to the direct long-ranged interactions 
between the two substrates, while Ai :S (T) > is the 
corresponding term generated by the long-ranged tails 
of the substrate potentials acting on the fluid particles. 
According to Appendix [A] and with J s,s (r — r') = J s /|r — 
r /|d+CT f or j-jjg interaction between substrate particles one 
has 



MT) 



<j{a-l) T(^±2 
and according to Appendix 



(2.16) 



AiAT,p) = 



~2v sPb (T,p)/(a-l) 
87r (d-D/2 r(l±£) 
a(a-l) r(^) 



(2.17) 

J l > s p b (T,p)p s (T)>0, 
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where (see Appendix [Aj 



err 



(2.18) 



Note that Ai < and ^4 S < 0, while Ai >a > 0, and that 
A s is an analytic function of T. If both the fluid and the 
substrates are governed by van der Waals forces, the con- 
stants Ai, Ai s and A s are not independent. As pointed 
out in Ref. [50] for such systems the parameters J 1 , J^ s , 
and J s can be expressed in terms of the corresponding 
atomic polarizabilities ai and a s , i.e., J 1 ~ af, J s ~ af., 
and J z ' s ~ aia s . (Note that this property implies that 
J /af — J s /af = J l ' s /aia s is independent of the par- 
ticular atomic or molecular species.) Therefore, one can 
write the sum of Ai, Ai iS , and A s as a perfect square, i.e., 



Ai + A, s + A S = - 



4n (d-l)/2 r( l±a) 



a(a-l) T(^) 
[J l p 2 b (T,fx) - 2J l > s p b (T,»)p g (T) + J S P 2 S (T)] 
~ -(a sPs - ai Pb ) 2 <0, (2.19) 

which leads to the conclusion, that for a Lennard-Jones 
fluid between two identical Lennard-Jones walls the ef- 
fective dispersion interaction for the film is always attrac- 
tive. Coating the substrate surfaces of the system with 
some additional material does not change this leading- 
order L-dependence and therefore does not change the 
above property. This result is in full agreement with the 
Dzyaloshinskii-Lifshitz-Pitaevskii theory [Hll H2] which 
also predicts that the effective interaction between iden- 
tical half-spaces separated by a thin film is always attrac- 
tive m [Hi . 

The only quantity in Eq. (|2.14p . which still has to 
be identified for our model, is the value of the coupling 
constants According to Refs. dSHEil 



I = v a /v 2 , 



(2.20) 



where v a and V2 are coefficients in the Fourier trans- 
form J(k) = ^ r exp(zk • r)J(r) of the interaction J(r) 

(see Eq.{2HD). One has 0, El J(k) = 
.7(0) [1 - v 2 k 2 + v a k a - w 4 fc 4 + 0(k 6 )}. It is easy to 
check that the short-ranged part </ sr (k) of the Fourier 
transform of the interaction has the form J sr (k) 

./;..; i + 2d - k 2 + o(fc 4 )] = j sr (o)[i - vfk 2 + o(k 4 )}, 

with J sr (0) = (1 + 2d)J l SJ and vf = 1/(1 + 2d) while for 
the Fourier transform of the long-ranged part Ji r (k) = 



Jlr(0)[l 



v%k 2 



v l £k 4 + 0(k 6 )] an analytical 



expression in closed form can be obtained only for the 
product J\ T (0)v^: 



J{0)v a = J h .(0)v l * = J l n{d,a), 



with 

n(d, a) 



d/2+l 



2" sin(7rcr/2)r[(d + a)/2]r(l + a/2) 



(2.21) 



>0, 
(2.22) 



for 2 < a < 4. This leads to 

A n(d, a) 



I = 



with 



1 + A j£(0)u£ 



J or 



(2.23) 



(2.24) 



while J/ r (0) = j\ T {Q)/J l is the ground state energy 
J h .(0) = J l Y,r 1/(1 + |r| d+CT ), where r with |r| = 0,1 



are omitted from the sum (see Eq. (|2.9p ). of a system 
with purely long-ranged interactions measured in units 
of J . Since both Ji r (0) and uJf depend also on the prop- 
erties of the interaction at short distances it is clear that 
for the precise determination of the values of I the use of 
numerical methods is unavoidable. It can be shown that 



4(o)t£ 



i 



2d 1 



\d+cr 



9(\r\ 



1) 



(2.25) 



For the "genuine" van der Waals interaction d = a = 
3 one has n(3,3) = 7r 2 /12 ~ 0.822 and if the system 
is discretized in terms of a simple cubic lattice one has 
J£(0)4 r - 1-692. In this case for A = 0, 1/2, 1, and 2 
one has I = 0, 0.222, 0.306, and 0.375, respectively. Note 
that for the Lennard-Jones (6,12) potential 



wlj = 4e 
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(2.26) 



the position of its minimum, which within our approach 
can reasonably be taken as the lattice spacing a, is at a = 
\[2 (To. The value — e of its minimum corresponds to J l sr . 
Thus the leading, attractive part — 4e(<Jo/r) 6 of the long- 
ranged interaction equals —2e(a/r) e , which leads to the 
conclusion that for pure Lennard-Jones potentials A = 2. 

We recall that for L ^> L C rit! u s > 0, and > 0, in 
line with Eq. I|2.14p the behavior of the Casimir force 
can be decomposed according to Eq. (|1.3p . The aim of 
the current study is to determine the leading finite-size 
behavior of the Casimir force in any region of the ther- 
modynamic parameter space of the model. As stated 
in Sec. HI for any y, X sl '(x,y) decays exponentially for 
|x| — * oo, while in that regime the other scaling func- 
tions X\ T and X\ r decay according to a power law. Based 
on the results of Refs. [13, IH, [23, HO] one expects that 
the deviation from the leading short-ranged behavior sets 
for L > £ln(£/£o)- For the scaling function X] r one ex- 
pects X\ T (x — > oo, 0) = Xix n ~ 2 , where X\ is a universal 
constant. 

Some information is also known about the asymp- 
totic behavior of X\ Y [x -> oo,0) (H^. In Ref. (HH| the 
authors have considered the behavior of the solvation 
force in a fluid with short-ranged forces (i.e., J 1 = 0) 
under the influence of long-ranged substrate potentials 
decaying at large distances z as Sv s z~° '. In the case 
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FIG. 3: (a) The typical behavior of the total Casimir force (per area of the film cross section, in units of feT, and without 
the direct substrate-substrate interaction) at A/i = as a function of the temperature scaling variable in a system with 
substrate-fluid and fluid-fluid long-ranged interactions. One observes that outside the critical region the force is repulsive and 
becomes attractive only near the critical point of the system due to the critical fluctuation of the fluid density. (If the direct 
substrate-substrate interaction is added to /total, the resulting force is attractive in the whole thermodynamic space.) The data 
correspond to L = 50 and s = A = 1 (where s = 1 for T = T c , see Eq. Q2.1ip ). The asymptote (see Eqs. I|2.15l) and l|2.17|t l 
is given by [3[a - 1) (A ttS + At) = 4sKp b + XK {47r d_1 r[(l + a)/2]}/{aT[(d + a)/2]}p b (1 - p b ) where in the present mean-field 
case d — a — 4, K = PJ 1 BI , and for K c see Eq. I|3.6p . Note that above the critical point p b = p c — 1/2, while for T < T c p b 
increases upon a decrease of the temperature. This leads to a clearly visible asymmetry in the asymptotic behavior of the force. 
(6) The same as in (a) but as function of the field scaling variable and at T c . One again observes that everywhere outside the 
critical region the total force (without the direct substrate-substrate interaction) is repulsive and becomes attractive only near 
the critical point of the system due to the critical fluctuation of the fluid density. 



Ai — A s = Ai iS = (i.e., for T > T c in a magnetic 
system in the absence of an external bulk magnetic field) 
they found numerically together with analytical mean- 
field arguments that the leading size dependence of the 
force / is of the order of L~( a+1 \ If valid this would 
lead to X*(x -> oo,0) = X s x d / 2 - 2+r >/ 2 . For a mean- 
field model with d = 4 and 77 = this would imply 
X 1 /(x — »• oo,0) = X s , i.e., this scaling function would 
tend to a constant away from the critical point (or vary 
~ lnar). Moreover, consistency with the structure of Eq. 
(|1.3p imposes the temperature dependence ~ t@~ u of the 
term oc L~( a+1 K This differs from the result in Ref. [H3 |. 
where the authors analyzed the system away from the 
scaling regime and purportedly reported that the corre- 
sponding term varies oc t^ 1 L~ l - a+1 \ 

Finally, we finish this section by recalling that if L < 
-^crit (or if u) a < or u)\ < 0) an expansion of the type 
given in Eq. (|1.3|l is not possible. In that case the effects 
of the van der Waals interactions will be important ev- 
erywhere in the critical region of the system, even at the 
bulk critical point. 



III. FINITE-SIZE BEHAVIOR OF FILM FREE 
ENERGIES AND OF THE CASIMIR FORCE 

In this section we investigate the finite-size behavior 
of the total Casimir force within the full temperature 



range. Experimentally this force is accessible as the so- 
called solvation force. In Subsec. IIII A I we focus on the 
critical behavior; the finite-size behavior off criticality is 
discussed in Subsec. IIII Bt Before passing to a detailed 
analysis of the influence of the different parameters on 
the behavior of the force, we first provide a view on the 
typical behavior of the total force (s = A = 1) as a func- 
tion of L/£t at A/i = (i.e., at bulk coexistence) (see Fig. 
G^a)), and as a function of L/^ at the critical temper- 
ature T — T c (see Fig. EJb)). The results are obtained 
within the model presented in the previous section choos- 
ing a film thickness L = 50 layers. (Compare Fig. [2] 
for L = 500; however, here the forces are not normal- 
ized by their value at T c .) One observes that the force 
(neglecting the contributions from the direct substrate- 
substrate interaction so that the conclusion stated after 
Eq. I|2.19p does not hold for /total considered here) is 
always repulsive outside the critical region while within 
the critical region it becomes attractive due to the critical 
fluctuations of the density. One observes in both cases 
for L/l- > 30 (where £ is either £t or £ M ) that the con- 
tribution {3 (a — l){Ai + Ai tS ) (indicated as "asymptote") 
provides a good approximation of the total force. Using 
Eqs. and ([2?L7j) . and p c = 1/2, it can be shown 

that 

p(a- l){Ai + A ltS ) = 4s Pb + — T [ ^ XKp b (l- Pb ), 

(3.1) 
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where K = (3J l SI . Note that, if s > and A > 0, (3(a - 
l)(Ai + Ai iS ) > 0, i.e., the force away from the critical 
region is repulsive. 

In the following subsection we analyze in detail the 
behavior of the force within the critical region and study 
its dependence on s, A, and L. 

A. Critical behavior 

Within the present mean-field approach, one has (3 = 
v = 1/2, i] = 0, and 7=1. Apart from logarithmic 
corrections, these mean-field values hold for d = d c = 4, 
provided the interactions governing the system are not 
too long-ranged. Since for genuine van der Waals inter- 
actions d = a, in the following we adopt d = a = 4, 
which leads to: 

f(t,H,L) ~ L- 4 {x[L/Z t ,(f3An)(L/Z 0tfl f, (3.2) 

WSo) -2 ,*^/^)" 1 ] 

+3/3 [Ai(T, n) + A ljS (T, fj,) + A S (T)} }. 

Note that for g? = ct = 4, 77 = the requirements a > 2— 77 
and a > (c?+2 — 77) /2 for the irrelevance of the long-ranged 
van der Waals type contribution to the behavior of the 
force are fulfilled [H(|. For the model considered here a 
direct estimate of L cr ;t yields L cr it = 32s (taking into 
account that here all distances are measured in units of 
the lattice spacing). Since s = 1 is the largest value of 
s for which we shall provide numerical results, we con- 
clude that for systems with L 3> 32 layers the van der 



Waals interaction is expected to give only corrections to 
scaling while for smaller L its effects will be important 
everywhere including the critical region of the system and 
will affect even the determination of the Casimir ampli- 
tude. Of course the above estimate does not tell the exact 
meaning of "much larger than" 32, but we expect that a 
factor of 10 should put one onto the safe side. Thus the 
expectations is that for L > 300 the van der Waals inter- 
action will provide only corrections to scaling within the 
critical region. Our subsequent analysis will demonstrate 
that this expectation is indeed valid. We shall consider 
films with thicknesses L = 50, 100, 500, and L = 1000 
layers. It will turn out that for L = 500 and L = 1000 
the van der Waals effects are small, whereas for L = 50 
and L = 100 it will turn out that they affect the be- 
havior of the force for all values of the thermodynamic 
parameters. 

In order to determine the scaling function X of the 
Casimir force one has to solve numerically Eq. I|2.6p for 
(j)*{r) and to insert this solution into Eq. <|2 . 5(1 . where 
p* = {1 + <j)*)/2. Due to the symmetry of the system 
one has 4>*{r) = 4>*{z). We determine the force from the 
lattice version of Eq. (fTTTj) (see Eqs. (|2T2|) and ([21410 : 

f(L,t\l,s = 0) = -|[ Wex (L + l,t|Z,» = 0) 

-w ex (£-M|Z,s = 0)] . (3.3) 
(There are alternative approximations of duj ex /dL which, 
however, yield the same leading behavior in L we are 
interested in.) 

According to Appendix [Cj for d = a = 4 the equation 
for the order parameter reads (see Eq. (|C15|0 : 



arctanh[(/>*(z)] = ^9 [Ap - AV(z)} + K J a 4 4>*{z) + a?" [<j>*{z + 1) + <f>*(z - 1)] + A ^ g 4 (\z - z'\)<t>*{z') 



z'=0 
|*'-*i>2 



where K = j3J l st , A = J l /J l SI (Eq. (|2T24j) l. and (see Eq. (fClHl ) 

g 4 (a) = tt 3/2 / dtt 3/2 E 4A (-t 4 )exp[-ta 2 ] 
Jo 

(V2a 2 
1 + f V2a 2 - l" 2 



7T 

23/4 



iV4 r 



1 

1/4 



sin ( iarccot[\/2a 2 + 1] ) — cos ( iarccot[\/2a 2 + 1] 



sin ( -arccot[\/2a 2 — 1] ) + cos ( -arccot[\/2a 2 — 1] 



r 



In Eq. (jaUJ a 4 = 7 + A(c 4 - 4), a\ n = 1 + A(cJ n - bulk system is 
1/2), where c 4 = 4 900 and c ff" = 1.028 are constants 
evaluated in Eqs. I|C17P and <|C 18|) . respectively. It is 

straightforward to show that the critical coupling of the K^ 1 {\) = 04(A) + 2a 4 m (A) + 2\^^g 4 {z), 

z>2 



(3.4) 



(3.5) 



(3.6) 
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so that the bulk critical point coordinates are (K — 
A' C (A), A/i = 0). Note that the position of the critical 
point depends on A. 

For (+, +) boundary conditions one has p(0) = p(L) = 
1, or equivalently 0(0) = <f>{L) = 1. We have imposed 
the boundary conditions in such a way that the number 
of liquid layers with independent degrees of freedom is 
L — 1 . (We recall our remark in the Introduction that the 
Casimir force hinges on the definition of what is called 
the film thickness.) The function 54(a) accounts for the 
fluid-fluid long-ranged interaction. We note that in many 
analyses 54(a) is typically neglected; even in numerical 
treatments of van der Waals systems the range of the 
interaction between the fluid particles is often actually 
truncated rendering effectively a short-ranged potential. 
One can cast Eq. (|3,4p also into a continuum form which 
is a close analogue of the one usually derived from a 
Ginzburg-Landau type Hamiltonian. The corresponding 
generalized Euler-Lagrange equation for the profile then 
reads (see Eq. (JC20j)): 



0*(z) + i(0*(z)) 3 = i/3(A M -A^)) 



+A^a 4 0*(z)+a™ 



20* (z) 



d 2 4>*(z) 



dz 2 



+ / gi (\z-z'\)<p*(z')dz' 
Jo 



(3.7) 



We start our evaluation of the Casimir force by first 
reproducing the analytically exactly available mean-field 
results for the force for (+, +) boundary conditions. That 
serves also as a valuable check for the model we have 
described above. 

For the following, Eq. l|3.4p defines the model the 
finite-size behavior of which will be investigated in detail. 
In this model the long-ranged substrate fluid interaction 
is present via the substrate potentials AV, while the long- 
ranged fluid-fluid interaction is reflected by A 7^ 0. Ac- 
cordingly, the fully short-ranged model corresponds then 
to AV = and A = 0. 

Before starting with the study of the influence of the 
different parts of the interaction on the behavior of the 
force we briefly comment on the universality of the scaling 
function X in Eq. l|3.2p within mean-field theory. 



1. Modified finite- size scaling for mean-field systems 

Here we assume that L is large enough so that one 
is allowed to consider only the scaling fields which are 
relevant in renormalization-group sense. In this case, 
within non-mean-field theories hyper-scaling and hyper- 
universality are valid so that the finite-size behavior of 
the singular part of the grand canonical potential density 
f3uj s = j3Cl s /V near the bulk critical point of the system 



(T = T c , fi = fx c ) is given by \5} 



l3u s (T,fi,L)=L- d X u 



L 



L 



6(A) '^(A) 



(3.8) 



and 



where £ t = ^(T - T±,/i = /i c ) = &(X)\t\ 
^ = ^(T = T c ,/i -» /i c ) = ^(X)\p c A^\-^ A are the 
second moment bulk correlation lengths. Here, as before, 
all lengths are measured in units of the lattice spacing a, 
t = (T — T c )/T c , Afj, = [i — fi c , and the parameter A 
reflects the dependence of the bulk system on the long- 
ranged component of the fluid-fluid interaction. Note 
that there is no nonuniversal normalization factor in front 
of the universal scaling function X UJ (xt,x fJ _). Indeed, 
due to hyper-universality lim T ^ T + /3uj s (T, fi c , oc)£f = Q, 
where Q is a universal constant; this is consistent with 
the limit L/£ t — > 00 of Eq. (|3.8j) with lim x _oc X u (x, 0) = 
Qx d . 

Within mean-field theory hyper-universality is lacking 
and generates a nonuniversal metric factor A(X), i.e., 



(3a Js (T,fi,L) = L-' i A(X)X i 



MF 



L 



6(A)' & (A) 



(3.9) 



It is easy to check that within our model the bulk grand 
canonical potential density /?w Si buik(7 1 , A 4 ) of the system 
is a universal function of /3/f3 c and (3Afi in the sense 
that /3ws,buik(r, /i) is not proportional to any A depen- 
dent term. Requiring this property to be compatible 
with Eq. 1|3.9|1 in the limit L/£t — * 00 , it follows that 
A(X) has to be of the form A(X) — Qmf [Co" (A)] 4 , where 
Qmf is universal (i.e., independent of A). Furthermore, 
one has X^f F (x — > 00, y) — x 4 X^ LF (x/y). Moreover, 
since /3w Sj t, u ik(r, fi) does not depend explicitly on A it 
follows that the ratio does not depend on A, which 

for our model can indeed be verified (see the following 
subsection). Thus Eq. I|3.9p turns into the form 



/3w s (T,/x,L) 



L 



X 



MF 



L 



L 



^{X)J Ut(A)'UA), 

(3.10) 

where Qmf has been incorporated into the scaling func- 
tion X^f F . X^ F is "universal", in the sense that it 
does not depend on A (neglecting the contributions due 
to the irrelevant scaling fields). Naturally, the ratio 
0uj 3 (T, fj,)/p c u] s (Tc, He) defines a scaling function that 
also does not depend explicitly on A (see Eq. Q3,9p ). 
but this ratio does not offer the possibility to discuss 
the dependence of [3ui s (T, /z) on other parameters of the 
model near the critical point because by construction at 
the critical point this ratio equals 1. 



2. The Casimir force in systems with short-ranged 
interactions. 

For (+, +) boundary conditions and within mean-field 
theory the analytical form of the scaling function of the 
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FIG. 4: (a) Comparison for A/i = between the behavior of the total Casimir force in our mean-field lattice model with short- 
ranged interaction (red dots) and the singular part of the force obtained within continuum mean-field Ginzburg-Landau theory 
[13] (black line). Both curves are normalized by their absolute values at the bulk critical point. This normalization is necessary 
because the Ginzburg-Landau theory predicts the force only up to a nonuniversal multiplicative factor. The lattice data have 
been calculated for L = 500. Both curves overlap perfectly, which demonstrates that the corrections due to the corrections to 
scaling, regular contributions, and the existing arbitrariness in the definition of L can be safely ignored if L is of order the of 
500 and the interaction is short-ranged. Note that the minimum occurs above the critical point at x t = sgn(t)L/^ t — 3.75. (6) 
The same as in (a) for T = T a as a function of the field scaling variable = the results for the singular part of the 

force within a continuum mean field Ginzburg-Landau type theory are from Ref. [21J]. The minimum of the force occurs at 
x M ~ —8.4, i.e., on the vapor side of the bulk coexistence curve (with the walls preferring the liquid phase). The force minimum 
is about a factor 10 deeper than the force at T c and ca. seven times deeper than the minimum along the temperature scan. 



Casimir force in systems with short-ranged forces has 
been obtained in Ref. flU ]. The corresponding result 
for Afi — is 

(i) X sr (y>0) = -[2^(fc)] 4 fc 2 (l-fc 2 ), (3.11a) 

with y = [2K{k)] 2 {2k 2 - 1), 

(ii) X sr (0>y>-7r 2 ) = -4^ 4 (A:), (3.11b) 

with y = [2K (k)] 2 (2k 2 - 1), 

(hi) X 3I {y<-Tr 2 ) = -4K\k){l-k 2 ) 2 , (3.11c) 

with y = — [2K(k)] 2 (k 2 + 1), where K(k) is the complete 
elliptic integral of the first kind, < k < 1. In Eq. (|3. 
the scaling variable of X SI (y) is y = tL x l v = tL 2 (here 
L is measured in units of a, i.e., in the scaling variable 
y L is dimensionless) and enters implicitly via y — y(k) 
which can be inverted uniquely to k — k(y) because y 
is a monotonic function of k. We note that X sr (y) is 
analytic for all values of y, because the film critical point 
(T C (L) , An c (L)) is located off coexistence at Afi c (L) ~ 
L- k / v ~ L~ 3 (see Fig.©) (MSI]. Obviously one has 
y > if A: > 1/ y/2 (with k = 1/ y/2 corresponding to the 
bulk critical point), > y > -tt 2 if l/%/2 > k > (with 
k = corresponding to the actual film critical point), 
and y < —it 2 if — 1 < k < (negative k describe the 
region below the bulk critical point). 

In order to check universal aspects of our model and 
the reliability of our numerical procedures, for A// = 
and T close to T c in Fig.[4ja) we compare the results for 



X sr (x t = sgn(t)Z/£ t , afp = sgn( A/x)L /g M ) = 0) obtained 
within the lattice model (Eq. IpTHI with Afi = AV = 
^rog = J 1 = 0) with the above analytical results of the 
continuum theory (l^ |. In addition, for T — T c and A/j, ^ 
in Fig.lUjb) we present the scaling function of the force 
X sr (x t = 0,Xfj,) as obtained within the lattice model in 
comparison with the numerical results of the continuum 
theory [2l[. Here £t is the bulk correlation length £(t — > 
±0,A^ = 0) = £o\t\- v , with v = 1/2, while ^(t = 
0,A^ -4 0) = Zo, M \An\- v/A with u/A = 1/3. For the 
continuum model £g~ = a, = a/v2 and £o,^ = a l 
[2ll . l58l | (where a has been introduced as a length scale in 
order to achieve compatibility with the lattice model). 

The Fourier transform g(k) of the density-density cor- 
relation function 



g(ri - r 2 ) 



has the form 



5(k) 



(p(r 1 )p(r 2 ))-(p(r 1 ))(p(r 2 )) (3.12) 

i^Cn)^^))-^))^))] 



(3.13) 



[p(l-p)]-i+0w(ky 



From Eq. (|3.12[) . using — 4J(k) = w(k) and the small 
k expansion of J(k), one obtains for the second moment 
correlation length £ (compare with Eq. (6.9) in Ref. (58| ) 



/'t'2 



4p(l-p) 



- 1 



T/T e 



(3.14) 
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FIG. 5: Comparison between the behavior of the total Casimir force normalized by its absolute value at the bulk critical point 
for a lattice model with short-ranged interactions for L = 50 (red dots) with the corresponding continuum Ginzburg-Landau 
result (black curve), (a) For A/i = there are small differences between the two curves for |sgn(t)L/£t| > 5 where corrections 
to scaling and the regular contributions start to show up. For \sgn(t)L/£t\ 3> 1 both curves decay exponentially and therefore 
they cannot be distinguished for jsgn(t)L/£t| > 10. (6) Same as (a) for T = T c and as function of Afj, in terms of the field 
scaling variable = L/^, The deviations are largest around the minimum. 




sgn(t) L/f t (A) x t =sgn(t) L/ft(A=0) 

FIG. 6: Comparison between the behavior of the total unnormalized (a) and normalized (6) Casimir force in a lattice model 
with short-ranged interactions for different values of L. While in (a) the curves for L = 500 and L = 1000 practically coincide, 
which means that the nonuniversal corrections are negligible, the corrections become apparent for the film thicknesses L = 50 
and L — 100. At the bulk critical point these corrections amount to 15% for L = 50 and 8% for L — 100. (b) If the curves 
from (a) are normalized by their values at the critical point they practically fall on top of each other and coincide with the 
continuum result. This occurs because the force is strongest at and near T c and, by construction, the curves for different L are 
forced to coincide at T = T c thus shifting the possible deviations from each other towards large values of the scaling variable 
for which the force decays exponentially so that in that range all curves again coincide with each other. 



Solving the bulk equation 4> = tanh[0(/3//3 c ) + /3A/i/2] 
for the order parameter one obtains from Eq. (|3.14|) that 
in the lattice system under consideration 

(3.15) 

For a system with a fluid-fluid interaction given by Eq. 
(|2.9p with J 1 = (i.e., short-ranged interaction) and dis- 
cretized on a cubic lattice one has V2 = 1/(1 + 2d), so 
that for d = 4 one has ^Jv^ ~ 1/3. Figures SJa) and 
SJb) demonstrate that the normalized scaling functions 
as obtained within the continuum model and within the 
lattice model practically coincide once the amplitudes 
£0 and £o,/j for the two models are chosen as described 



above. Thus on the mean-field level the present lattice 
model is indeed suitable for investigating the Casimir 
force. We emphasize that at T = T c the minimum of 
the force does not occur at bulk coexistence /1 = /i c 
but at a nonzero value A/i, determined for each L via 
= sgn(A/x)i/^ jU ~ —8.4, i.e., on the gas side of the 
bulk coexistence curve (with the walls preferring the liq- 
uid phase) . Similarly, if A/z = the minimum of the force 
is not at but above T c , given by xt = sgn(t)L/^ t — 3.75. 
We note that Fig. [4] actually provides a comparison be- 
tween the behavior of the total Casimir force in a lat- 
tice model with short-ranged interaction with the singu- 
lar part of the force obtained within continuum mean- 
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FIG. 7: Comparison between the behavior of the critical Casimir force normalized according to Eq. (13. 19[) for different values of 
L, A, and s. The strength of the long-ranged fluid-fluid interaction is reflected by the value of A (the larger A the stronger is the 
long-ranged tail of the interaction in comparison with the short-ranged part (Eq. (J2T24J) ) ) and the strength of the long-ranged 
substrate-fluid interaction is reflected by the value of s (the larger s the stronger is the effective interaction with the substrate 
(Eqs. I|2.8p and l|2.1ip )). We consider systems with L = 1000 (a), 500 (&),100 (c), and L — 50 (d) layers. If L decreases one 
observes a more pronounced dependence on s and A. While for L = 1000 and L = 500 all curves are very close to each other, 
they significantly deviate from each other for L — 100 and L = 50. For all L the maximal deviations with respect to the 
short-ranged case corresponding to s = A = occur for s = 1, A = (which strengthens the force) and s = 0, A = 1 (which 
weakens the force). As nonuniversal effects they depend on the absolute value of L. While for L = 1000 the maximal deviation 
is (still) about 10%, for L — 50 it reaches 100%. This profound dependence on the values of s and A is most visible for L — 50. 
Contrary to the situation for L = 1000 and L = 500, for L = 50 none of the curves are even close to each other. The data show 
that the dependence on s is more sensitive than that on A: changing the value of s leads to more pronounced deviations from 
the results for short-ranged forces than changing the value of A. Furthermore, for smaller L the dependence on s is stronger 
and the curves vanish less rapidly for x t — > oo. We recall the asymmetry in the asymptotes of the total force for x t — > — oo and 
xt — » oo (see Fig. [3ja)) . In order to obtain its critical contribution the asymptote for T > T c was subtracted. For this reason 
one observes for /^n T a tendency to increase for x t — > — oo. 



field Ginzburg-Landau theory. If L is small enough the 
corrections to the universal behavior, which are due to 
the corrections to scaling, to regular contributions, and 
to the existing arbitrariness in the definition of L (see 
below as well as Appendix B in Ref. [!§]), will become 
visible. For L = 50 this is shown in Fig. [5l The role 
of L is visualized in Fig. El Since the Casimir force is 
a derivative of the excess free energy with respect to L 
it depends on the definition of L for a given system. In 
the current analysis this distance is taken to be the lat- 
tice spacing times the number of layers with independent 
degrees of freedom (which is L — 1) plus two times half 
the distances between the outermost layers with frozen 
degrees of freedom and their adjacent inner layers, i.e., 



La. 



3. The Casimir force in systems with van der Waals type 
interactions 

In this subsection we analyze the effect of the range of 
the fluid-fluid and the substrate-fluid interactions on the 
behavior of the Casimir force, i.e., we consider systems 
in which either A or s, or both, are nonzero. Note that 
A ^ implies I ^ (see Eq. p33]l ). To this end we first 
determine the bulk correlation length amplitudes £o and 
£o,m which depend on A. Due to Eq. I|3.15p this requires 
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to calculate the Fourier coefficient V2, for which we obtain Figs. [7] and [9] show the temperature and field dependence 



V2 



_ M 2 

2d Z^r 1+jrp 



E 



l + 2d + A£ rTT ^#(|r|-l) 



(3.16) 



For d = a = 4 and for a hypercubic lattice a numerical 
evaluation yields 



''2 



1 + 0.829 A 
9 + 2.152 A' 



(3.17) 



For A = 0, 1/2, 1, and 2 one thus obtains £q = y/v2 = 
1/3, 0.375, 0.405, and 0.447, respectively. 

The corresponding numerical results for the behavior 
of the critical Casimir force / cr it are summarized in Figs. 
[3 El El and[l0l (/crit is that part of the total force /total 
calculated with the regular background term subtracted; 
this procedure corresponds to the analysis of experimen- 
tal results if one subtracts from the value of the force 
measured around T c the asymptote obtained by measur- 
ing this force well above T c .) While Figs. [7] and [8] deal 
with the temperature dependence of the critical Casimir 
force at coexistence, Figs. [9] and [10] visualize the depen- 
dence on the excess chemical potential at T = T c . In 
Figs. [8] and [10] the force is normalized by its value at 
the bulk critical point. This latter presentation does not 
allow one to compare the forces for different s and A at 
the critical point, because there by construction this ratio 
equals 1 for all parameters. This partial lack of insight 
from this ratio can be overcome by choosing an alterna- 
tive normalization, which keeps the dependence on A and 
s even at the critical point: 

/crit (x t , x^ | A, s) (3.18) 

- /crit (it, ^|A, s)/|/crit(0, 0|A = 0, S = 0)|. 

Here / C rit(0,0|A = 0, s = 0) is the leading asymptotic 
behavior (d — 1) A +:+ L _d of the critical Casimir force in 
a system with short-ranged forces. This can be inferred 
from experimental data of actual systems with dispersion 
forces by considering the limit L » L cr it (Eq. I|1.7p ) 
of the difference between the total Casimir force /total 
and the regular background contribution (see Eq. (|3.ip 
and the asymptote in Fig. [3]) at the bulk critical point 
(T = T c ,A/i = 0). In theoretical analyses / cr it(0,0|A = 
0, s = 0) can be determined directly by switching off the 
long-ranged parts of the interaction, i.e., by taking A = 
and s = from the outset. As explained in Subsec. 
MI All within mean-field theory the rhs of Eq. (f37T8j) 
is not a ratio of universal scaling functions (Eq. I|3.9p ) 
as it is the case for the rhs of Eq. (|3- 18|) beyond mean- 
field theory, i.e., for d < 4. However, on the basis of the 
constructions in Subsec. MI A H one can find a mean-field 
estimate for /critC^tj %n\\ s) in d — 3 by determining, 
e.g., within our present mean-field approach, 

/™r(^*,iA, S ) 

= {[/crit(xt,a; A1 |A, S )/|/crit(0,0|A = 0, S = 0)|] x 
[eo + (A = 0)/£+(A)] 4 } MFT (3.19) 



of /™f T , respectively, for L = 1000, 500, 100, and 50. 

For decreasing values of L one observes that the de- 
pendence on s and A becomes more pronounced. The 
maximum value of s we have considered is s = 1. As 
stated in the beginning of Subsect. MI A| in this case one 
expects that for L > 300 the van der Waals interactions 
give rise only to corrections to the leading short-ranged 
behavior and that for thinner films there are significant 
deviations in the whole range of values of the scaling vari- 
ables. These expectations are confirmed by Figs. [7] and 
[9l While for L = 1000 and L — 500 all curves are very 
close to each other they deviate from each other signif- 
icantly for L = 100 and L = 50. Note that in all cases 
the force has a minimum above the critical point. The 
minimum is deepest for s = 1, A = and the force is 
weakest for s = and A = 1. 

Figures [8] and [10] show the same mean-field data 
as in Figs. [7] and [H respectively, but now normal- 
ized by the actual values at the bulk critical point, 
i-e., fcTit(xt,x^\X, s)/|/ cr it(0,0|A, s)\. Whereas the coin- 
cidence of the curves at Xt = and x^ = 0, respectively, 
is enforced by construction, for L = 50 and L = 100 
there are significant differences for \x t \ > 10, but barely 
affecting the minimum at x t — 3.75 (Fig. [8]), and for 



< —4 including the minimum at x^ 



-8.4, but 



barely affecting the scaling regime x^ > 0. The occur- 
rence of the force minimum at negative values of x M , i.e., 
on the gas side of the bulk gas-liquid coexistence curve, 
can be understood by recalling that the Casimir force is 
a fluctuation induced force. The fluctuations are max- 
imal at a small negative bulk field which is needed to 
neutralize the ordering effect (which suppresses the fluc- 
tuations) due to the substrate-fluid potential and thus of 
the (+, +) boundary conditions. The ordering effect of 
the surface potentials can also be neutralized by increas- 
ing the temperature to values slightly above T c . There- 
fore the magnitude of the force attains its maximum for 
T = T c at a negative value of x^ and for A/j, = at a 
positive value of Xt . 



B. Finite-size behavior off criticality: capillary 
condensation regime 

To a certain extent Fig. [3] has already discussed some 
features of the force outside the critical region, where the 
total force is repulsive. Off criticality the force is well ap- 
proximated by (3{(J — l)(Ai + Ai tS ) given by Eq. (|3.1[) . 
The small deviations A/ of the force from its asymptotes 
(see Fig. ([3])) can be explained in terms of an effective 
finite-size contribution A/j,l to the excess chemical po- 
tential A/i, which scales as 1/L, i.e., the confined fluid 
at A/j, has approximately the properties of a bulk fluid 
at Afi + Afi L [5i|. This implies (p = (1 + <j>)/2) that 
(f> ~ (j) b {A^ + A/ij,) ~ (f> b {Afi) + Xh(A y u)A / u L , with the 
second term producing a contribution A/ ~ 1/L rela- 
tive to the leading first term. We recall that the total 
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FIG. 8: Same data as in Fig. [7] but normalized by the value of the critical Casimir force at the bulk critical point. This 
normalization enforces coincidence of the curves at Xt = 0. The figure demonstrates that this kind of normalization suppresses 
the visibility of the influence of the long-ranged interactions, which is actually rather strong for L — 50 and L = 100 (see 
Fig. [7)). Nevertheless, for such values of L and large scaling arguments the data demonstrate the dominant character of the 
finite-size contributions to the force which are due to the van der Waals fluid-fluid (A) and substrate-fluid interactions (s). 



force discussed in Fig. [3] has turned out to be repulsive 
outside the critical region because it does not contain 
the contribution A s from the direct substrate-substrate 
interaction. A s carries only a weak and smooth temper- 
ature dependence via p s (T) , which is not included in the 
degrees of freedom considered by the functional (|2.2p . If 
the contribution A s is added to /total, the resulting force 
is attractive (see Eq. (|2.19p ) throughout the whole ther- 
modynamics space. Besides the critical regime there is, 
however, one additional region where the approximation 
of the force by (3(a — l)(Ai + Ai ;S ) is also not valid and 
where the force can be attractive, too, even if A s is ne- 
glected. 

This is the so-called capillary condensation regime, 
which will be discussed in the following. This capillary 
condensation regime occurs at low temperatures T < T c 
and for LA/i = O(l) with Au < if the walls prefer 
the liquid phase [U, [6(| EH, [g^]. In accordance with 
Fig. Q] capillary condensation means that upon increas- 
ing the chemical potential in a film of thickness L there 
is a first-order phase transition between spatially inho- 
mogeneous gas-like and liquid-like configurations at val- 
ues A/i cap < 0, i.e., before reaching the bulk gas-liquid 
coexistence curve A/i = 0. The significant and interest- 



ing features of the force induced by capillary condensa- 
tion are displayed in Figs.Ql] and [321 For a film with 
(+,+) boundary conditions, for A/i < with |A/x| suffi- 
ciently small it is favorable for the system that (fiAfi < 0, 
i.e., the equilibrium phase to be liquid- like. This means 
that cfi is positive and thus follows the boundary con- 
ditions provided by the effective substrate potentials. 
If A// becomes more negative, at a given undersatura- 
tion A[i = A/x cap (T, L) the sign of the order parameter 
changes abruptly to a negative value following the bulk 
field and forcing the appearance of two interfaces within 
the system. Accordingly, for large L the corresponding 
contribution in SI which is linearly proportional to A/i 
changes by 20bA/iL, thus leading to a L-independent 
change 20{,A/i in the force. (This simple argument as- 
sumes that within this regime (A/z < 0, low T) in a large 
part of the film the order parameter profile <p(z) can be 
approximated by its bulk value (fib.) Since the jump in 
the sign of <fi occurs if the advantage of having (fib with the 
same sign as A//, i.e., </>&A/i > so that (fib < 0, compen- 
sates the disadvantage of creating two surface free energy 
contributions in the system with a negative order param- 
eter occurring inspite of (+,+) boundary conditions, it is 
easy to understand that A/i cap ~ 1/L. In Fig. [12] the 
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FIG. 9: Same as Fig. [7] as function of the field scaling variable a; M = sgn(A^t)L/£ M (A) for T = T c . The minimum of the force 
occurs for negative A/i, i.e., on the gas side of the bulk gas-liquid coexistence curve. The deepest minimum of the force occurs 
for s = 1, A = 0. For L — 1000 the force at the minimum is about 10 times stronger than the force at the critical point. If 
L decreases this minimum becomes even deeper; for L = 50 the corresponding force is about 15 times stronger than at the 
critical point. The force is weakest for s = 0, A = 1. Similar to Figs. [7]and[8]one again finds that the dependence on s is more 
sensitive than that on A: changing the value of s leads to more pronounced deviations from the results for short-ranged forces 
than changing the value of A. We recall the asymmetry in the asymptotes of the total force for x t — > — oo and xt — > oo (see 
Fig. [3](b)). In order to obtain its critical contribution the asymptote for T > T c was subtracted. 



total force is shown as function of the chemical potential 
A/z, For A[i > AjUcap, up to a leading order, it does 
not depend on L. As demonstrated numerically for large 
L only the maximum absolute value of the force is L- 
dependent, decreasing as ~ 1/L due to the L-dependence 
of A/x cap (T, L) . Figures [TT] and [32] present results for a 
system with short-ranged interactions (s = A = 0) and a 
system in which both the fluid-fluid and substrate-fluid 
interactions are long-ranged (s = A = 1). It turns out 
that the behavior of the force in the considered (T, fx) 
region depends on s and A mainly via the dependence of 
A/Xcap on s and A. This dependence is studied systemati- 
cally in Fig. [131 It shows that A^ cap depends linearly on 
K, L^ 1 , s and A. Furthermore the maximal amplitude of 
the force also increases linearly if any of these parameters 
increases. 



IV. SUMMARY AND CONCLUDING 
REMARKS 

In view of future experiments exploring the critical 
Casimir forces in classical one- or two-component flu- 
ids confined by parallel substrates with the same pref- 
erence for their coexisting fluid phases we have analyzed 
the actual and generic interplay between the long-ranged 
fluctuation induced force and the long-ranged dispersion 
forces acting both between fluid molecules and between 
fluid and substrate molecules. 

From general scaling arguments and explicit mean-field 
model calculations we have obtained the following main 
results: 

(1) Based on general renormalization group arguments 
(Subsec. II Bp we have established a relevance- 
irrelevance criterion (Subsec. II C|l for the impor- 
tance of dispersion forces within the critical region 
of the system. It turns out that if the thickness 
L of the liquid film is much larger that L cr i t (see 
Eq. (|1.7p ). within the critical region of the system 
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FIG. 10: Same data as in Fig. [9] but normalized by the value of the critical Casimir force at the bulk critical point. Whereas 
in Fig. [8] as function of xt the deviations occur for xt > 10, here the deviations are most significant for < —4, affecting in 
particular the minimum, but barely the scaling regime x M > 0. This is due to the fact that here the minimum of the force is 
rather deep (approximately 7 times deeper than the corresponding one in Fig. [8j so that the spread in the curves for x^ > 
is difficult to see, even for L = 50, on the scale used in the plot. 



the dispersion forces provide only corrections to the 
leading behavior of the force. However, outside the 
critical region, i.e., for temperatures T ^ T c or for 
undersaturations A/i ^ 0, the influence of the dis- 
persion forces becomes important. This is of exper- 
imental importance because it is difficult to thermo- 
dynamically position the system right at the criti- 
cal point. In the opposite case, i.e., for L < L ait 
the contribution of the van der Waals (dispersion) 
forces is always important and cannot be neglected 
even within the critical region of the system, in- 
cluding the bulk critical point. Within a mean- 
field type model (Eqs. pi!]) . flU}, and 
(|2.10p ) we have studied the behavior of the force 
as a function of the strength of the long-ranged 
part of the fluid-fluid interaction, reflected by the 
parameter A (Eq. l|2.24p ). on the strength of the 
long-ranged part of the effective substrate-fluid in- 
teraction, reflected by the parameter s (Eq. (|2.11|) ). 
on the thickness L of the film, and on the (relevant) 
temperature and chemical potential scaling fields. 
The interactions are taken to decay with the dis- 
tance r between the fluid particles as r^^^ and 
as z~ as for the distance z of a fluid particle from a 



single substrate. (Different from other theoretical 
model calculations we have taken the long-ranged 
tail of the fluid-fluid interaction fully into account, 
which amounts to respecting the fact that the func- 
tion gd{a) (Eq. (|3.5p and Fig. [14]) is nonzero.) For a 
d-dimensional system we consider the case a > 2— T) 
and a s > (d + 2 — rj)/2, which guarantees that the 
long-ranged tails of the interactions change neither 
the bulk nor the surface universality classes (3ll. I32T] . 
The modified finite-size scaling, which is needed if 
one - as done here - compares within mean-field 
theory systems with different A, is presented in Eq. 
(|3.10p . The results for the behavior of the force 
as a function of L/£t and L/^, where £t and £^ 
are the second-moment bulk correlation lengths for 
A/i = and at T c , respectively, are summarized in 
Figures [3, El El and [TO] These explicit results con- 
firm the relevance-irrelevance criterion mentioned 
above. The numerical results presented are for films 
with thicknesses L = 1000, L = 500, L = 100, and 
L = 50 layers in our simple cubic lattice model. 
We have varied A and s between (short-ranged 
model) and 1. Figures [3 [U [9] and fTOl show that 
for L = 1000 and L — 500 one observes universal 
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PIG. 11: Grand potential density in excess of the bulk con- 
tribution as function of L /3A/i for the two competing profile 
shapes (insets) in the capillary condensation regime for (+, +) 
boundary conditions; K = /3Jj r - For s — A = and s — A = 1 
the curves intersect at L/3A/i cap ~ —0.42 and —0.95, respec- 
tively. The thermodynamically stable states correspond to 
the minimum branches of /3u>. They are presented by the 
thick parts of the lines. The thin parts of the lines indicate 
metastable states. The full lines correspond to gas-like states, 
the short-dashed lines to liquid-like states. 



scaling, i.e., independence on the values of s and A, 
whereas this is lost for L — 100 and L — 50, which 
is in line with L < i C rit- The comparison between 
Figs. [7] and [Hon one hand and Figs. l9l and fTOl on 
the other hand shows which kind of analysis of the 
effects enhances or suppresses the visibility of the 
long-ranged forces. It turns out that L cr i t is sur- 
prisingly large (see Eq. (|1.7p ) which enhances the 
importance of taking into account the dispersion 
forces for the interpretation of experimental data 
of critical Casimir forces in fluid systems. 

(2) If all long-ranged interactions are taken into ac- 
count - the fluid-fluid, the substrate-fluid, and 
the substrate-substrate one, for the case of (+, +) 
boundary conditions studied here we confirm ex- 
plicitly within our model that the effective interac- 
tion between the substrates is attractive through- 
out the whole temperature range (see Eq. l|2.19p ). 
The latter is described by three different Hamaker 
constants which account for the aforementioned in- 
teractions. If the direct substrate-substrate inter- 
action is omitted (which amount to subtracting the 
force from the corresponding null experiment) the 
force is repulsive (see Figs. [2] and [3]) except within 
the critical and the capillary condensation regime. 

(3) For short-ranged systems, within the critical re- 
gion our lattice model reassuringly reproduces the 
universal scaling function of the critical Casimir 
force as obtained from continuum field theory, pro- 
vided L is sufficiently large (Subsec. IIII A 21 and 




-0.02 i ■ ■ ■ 1 

-0.02 -0.015 -0.01 -0.005 

FIG. 12: Total force (obtained from Eq. |T2~T2|l ) per unit 
area and in units of (3 as function of the chemical potential 
A/i < for /3\Afj,\L = 0(1) and T < T c . The discontinuity of 
the force due to capillary condensation occurs at A/i cap (T, L). 
For A/i > A/i cap (T, L), with A/i cap ~ 1/L, the leading be- 
havior of the force is basically independent of L and is given 
by 2(f>bAjj,. Note that the van der Waals contribution of the 
interaction strongly increases the maximal absolute value of 
the force. This occurs because both the long-ranged tails of 
the fluid-fluid (A) as well as of the substrate-fluid (s) inter- 
actions promote the ordered state and increase the absolute 
value of A/i cap beyond which the coexistence of a reservoir 
of gas with a liquid-like film is no longer possible. In this 
plot K = 0.2 corresponds to T/T C (A) = 0.556 and 0.453 for 
A = and A = 1, respectively. Both for A/i < A/i cap (T, L) 
and A/i > the leading behavior of the force is determined 
by the Hamaker term f3(a - l)(Ai + Ai tS ) given by Eq. I|3.ip . 
In these regions the magnitude of the force is of the order of 
L~ a (i.e., L~ 4 within our model), while in the capillary con- 
densed regime the maximal magnitude of the force is of the 
order of L~ x . For this reason on the present scales the curves 
for A/i < A/x cap (T, L) are not distinguishable. 



Fig. 2]). Surprisingly, even for short ranged forces 
(A = ,s = 0) L = 50 and L = 100 turn out to be too 
small in order to reach the universal scaling func- 
tion given by the continuum field theory; L = 500 
and L = 1000 are large enough in this respect (see 
Figs. [5] and EH). 

(4) At low temperatures, A/i < 0, and LA/i = 0(1) 
the force undergoes a jump as function of the excess 
chemical potential at a certain value A/i cap < 
caused by capillary condensation (see Figs. [H [TTJ 
and IT2]) . The position of the jump (Fig. |T2j) and 
its magnitude (Fig. [TJJ) vary linearly as function of 
1/T, 1/L, s, and A (Fig. [TJ). 
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FIG. 13: The curves illustrate the loci A/i cap (T, L\\, s) for 
capillary condensation leading to a jump of the Casimir force 
as a function of the coupling K — Js r /(fcsT) (a), the inverse 
linear size L (6), the relative strength of the surface-fluid po- 
tential s (c), and the strength of the van der Waals tail of the 
fluid-fluid interaction A (d). All these dependences turn out 
to be linearly decreasing within the studied parameter range. 
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APPENDIX A: DERIVATION OF THE 
EFFECTIVE HAMILTONIAN FOR A VAN DER 
WAALS SYSTEM CONFINED BETWEEN TWO 
PARALLEL PLATES 



Pj € {0, 1} and p\ € {0, 1} denote the occupation num- 
bers for the liquid and substrate particles, respectively. 
Since only the part {pj} corresponding to the fluid de- 
grees of freedom exhibit criticality at T c and the fluctua- 
tions of the substrate degrees of freedom {pj} are unim- 
portant in this respect, one can replace the latter ones by 
their mean-field values. If the fluid is in contact with a 
particle reservoir at a given (excess) chemical potential p, 
and temperature T, the partition function for the liquid 
is 



-/?U-/x£pJ 

4 4 E 
+2 ]T JlM + ^A 

i,jec jec I 



= E CX P 



(A2) 



where p\ = (pf). We assume that the solid substrate is 
only weakly influenced by the presence of its surface, so 
that pf = p s = const in S, which leads to 



Z = 



E ex P 



. jeC \ ies I 



+2/3 ]T J\jp\p\ 



(A3) 



We consider a fluid composed of particles interacting 
via a pair potential w l (r—r') = — 4j'(r— r') and bounded 
by substrates whose particles interact with the liquid par- 
ticles with a pair potential w l ' s (r — r') = — 4J^ s (r — r'). 
Within the lattice gas model for any given configuration 
C of particles {pf,p l j}, i £ S, j € C, where C and S 
denote the spatial region in a simple cubic lattice with 
lattice constant a occupied by the liquid and substrate 
particles, respectively, the energy of the fluid is given by 

E "•!:>:/'' • \ E <M 
= - 4 E - E < A| ) 



Modeling the pair potentials as 

4 tj = JilSQTi-rjD+SQn-Tjl-l)} 

+jl/(l + \r i -r j \ d +°) fl(| ri - rj |-l),(A4) 

(note that 0(0) = and 8(x) = 1 for x = and zero 
otherwise) and 

4? ee J^J(| ri - rj |-l) 

+J l ' s /\n - r,f+^ 0(\ Ti -Vj\- 1) (A5) 

one finds for J2i '■ 
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r 2l( rf +^)/ 2 



oo [(z, + n) 2 +r | + r | + . . . + r 2](<*+-*)/2 



00 f 00 27r( d_1 )/ 2 
tin / dr — 
Jo ] 

l)/2 r (~^) 



dri(zj + n) 



-<r.— 1 



r ( 1+2^ 

-D/2 1 \ 2 J 



' d+<r s \ "J 



(A6) 



where Zj > 1 characterizes the distance of the particle 
Pj from the boundary. We consider the fluid particles to 
be in the region < z < L, where L is the width of the 
film confined between the two surfaces. Therefore, the 
partition function is 



Z = E ex P 



2 E J ljPiVi 
i,jec 



, (A7) 



i.e., the system is equivalent to one with a spatially vary- 
ing chemical potential pj = p — Vj acting on a particle pj 
at a distance Zj + 1, < Zj < L, from the left boundary 
surface and at a distance (L + 1 — Zj) from the right one 
where Vj = V(zj) is the superposition 

V(z) = - Ps J l s f[S(z) + d(L-z)} 

+v s [{z + 1)-*° + {L + 1- z)-°°} , (A8) 



where 



•l+o-s 



PsJ 



l.s 



(A9) 



In the present study we choose J s ;. s such that p(0) = 
p(L) = 1, where p(z) — (p\). This is implemented by 



taking J s 



Ls 



oo, which is known as (+, +) boundary 



conditions applied to the system under consideration. 

Note that in Eqs. ((A6|l and ((A8|l only the leading z- 
dependence of the substrate potentials has been retained. 
We remark that there is some arbitrariness in what one 
calls the distance z. One might measure it from the last 
substrate layer or, say, from the midplane between the 
last substrate layer and the first liquid layer. These dif- 
ferent choices as well as taking into account the actual 
discreteness of the substrate lattice in normal direction 
(see the second line in Eq. i|A6|0 lead to additional terms 
decaying as z~ <Ts ~ 1 , z~ CTs ~ 2 , etc. In the following we do 
not consider such subdominant contributions. 

One can easily determine the (leading) L-dependent 

contribution Atoix^ in the excess grand canonical po- 
tential uj ex which is generated by the direct interaction 



of the two substrate half-spaces a distance L apart. Ac- 
cording to Eq. I|A8|) . one has 



E PsV s Z~ 
z=L + l 

_ 47r (d-l)/2. 



(A10) 



A S (T)L 



<r s r 



1 2 \ pU 

d-\-a s \ r & 



oo 

2 7S 1 dz; 



L+l 



where A s is given in Eq. I|2.16p . Also in this expression 
we have retained only the leading L-dependent part. 



APPENDIX B: FINITE-SIZE BEHAVIOR OFF 
CRITIC ALITY 



1. Finite-size contribution due to the long-ranged 
part of the fluid-fluid van der Waals interaction 



Away from T c the replacement of the fluctuating par- 
ticle density by its mean value, which is inherent to the 
present mean-field approach, is a reliable approximation. 
This leads to the following finite-size contribution Aa4x 
to the excess grand canonical potential per area A that is 
due to the long-ranged tail of the fluid-fluid interaction 
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A^oo 2A 



E 



(r lh z),(r' r z')eC 

-2J l pl £ 



w(r\\, z,r^ z')p(r h z)p(r' ll , z') 
' (Bl) 



' 1 + [r| + (z - z')' 1 



(d+a)/2 

r d ~ 2 dr 



MT,p)L-° +1 



z^z' 
l)/2 1 \ 2 J 



r(^) ji 



dz I dz' \z- z'\- (a+1) 



where A\ is defined in Eq. I|2.15p and the above inte- 
gration over z has been performed under the restriction 
\z — z'\ > 1, which mimics the regularization of w at 
r = r'. In all above steps in Eq. <|B1|) only the leading 
L-dependent terms, which are not part of the bulk or the 
surface contributions to the grand canonical potential, 
have been retained. 



2. Finite-size contributions of the long-ranged tails 
of the substrate potentials for T <C T c 



Away from T c the leading order finite-size contribution 
to the grand canonical film potential due to the long- 
ranged tail of the substrate potential V{z) = kJ],i s 5(z — 
1) + v s z~ a of, say, the left wall is 



Au,C«.«) c p b v s (j2 z ^ 



p b v s I dzz-°° = -j^-L-°° +1 , 



where we have again retained only the leading L- 
dependent part. 

Taking into account the contribution of the second 
surface, which renders a factor of 2 for the leading in- 
dependent part, we conclude that away from T c the con- 
tribution of the tails of the substrate potentials to the 
behavior of the Casimir force is 



A/ 



(l,s) _ 



dAtoj 1 ^ 
dL 



= -2(3p h v s L-°° = (3{a s - l)A 



La, 



r* (B3) 

which is basically the result derived in Ref. [5Q( (see Eq. 
(6.11) therein) and is in full accordance with the findings 
of Ref. [H3|. The constant A l:S is defined in Eq. (f2A7|l . 



APPENDIX C: ORDER PARAMETER PROFILE 
IN A FILM 



Since in the present study the external potential 
V{t) = V{r\\,z) = V(z) depends only on the coordinate 
normal to the confining plates and there is no symme- 
try breaking in the lateral directions, the density profile 
also depends only on z so that Eq. (|2.6p turns into the 
following equation for the order-parameter profile: 



4>*(z) — tanh 



+-f3(Ap J -AV(z)) 



(CI) 



where 



Jl(z) ee £ J(r,| -rj | ,,) = X: J(r h z). (C2) 



Measuring distances in terms of the lattice spacing and 
taking the interaction of the form (see Eq. I|2.9p ) 



J(r) = ^[«J(|r|) + J(|r|-l)] + 



J 



1+ r 



d+u 



M-i) 



(J> i -J')5(\r\) + (J^-~J l )5(\r\-l) 



rl jh 
J 1 

'l + \ T \d+a 



(C3) 



one can further simplify the sum on the right-hand side 
of Eq. JEU. Using the identity [63| (see also Ref. Q) 



1 f°° 

—— = / dt exp(-yt) e- 1 E at<x {-t a ), (C4) 

i + y a Jo 



(B2) where 



k=0 



a > 0, (C5) 



are the Mittag-Leffler functions, the sum over rn in Eq. 
(|C2[) can be rewritten as 

Ji(z) = {(J l sl .-J l ) + 2(d-l)U r -±J l )]6(z) (C6) 
+(•& ~ \J 1 )W - 1) + 5(z + 1)] + J^z) 

with 

P OC 

J l dA z ) = Jl / dt t^' 2 - 1 ^(-t^/ 2 ) x 

Jo ~ 

^e- tr ii | > (C7i 
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TTN^" 1 )/ 2 d+o 



The main advantage of the above form is that it factorizes reflects the contribution of the long-ranged van der Waals 
the summation over the components of ry. Using the tails of the interaction, and the function 
Poisson identity 

oo / — oo 

E ^ = S E e ^ 2,i2/t ( C8 ) 

V 

one has 



2 2/ . .2 

e~* n ft - ta (C12) 



Y j J l ^(z-z')4>(z') = J l L 

z'=0 I 



(C9) 



+ jr[g d A\z-A)+9rA\z-A)}^') 



= 



where 



^ 1 + r 



d+cr 

II 



(CIO) 



is a constant depending on d and cr. The function 



9dA a ) = I dt (^) ' -' 1 



/o 



7T\( d - 1 )/ 2 d+a , 



concerns mainly the coupling between neighboring sites 
in adjacent layers (see below). Here n 6 is a (<i— 1)- 
dimensional vector with integer components because all 
lengths are measured in units of the lattice spacing. It is 
easy to show that max t exp(— TT 2 n 2 /t + t(z — z 1 ) 2 ) is at- 
tained at t = 7r|n|/|z— z'\ and equals exp(— 27r|n| \z — z'\). 
Due to this exponential decay in the last term of Eq. 
(|C9P one is able to take into account only the terms with 
|n| = 1 and \z— z'\ = 1, which amounts to the approxima- 
tion g^'a( a ) — ffdCT^l) = ^da- ^ is straightforward to 
take into account additional, smaller terms corresponding 
to |n| = 2, 3, • • • and \z — z'\ — 2, 3, • • • , but it turns out 
that already the contributions stemming from |n| = 1 
and j 2 — z'\ = 1 are numerically very small. Thus the 
size dependent contributions due to the last term on the 
right-hand side of Eq. 1|C9|1 are exponentially small in L. 



J 



For d = a = 4 the corresponding Mittag-Leffler function can be expressed as 

1 



£4,4 M 4 ) 



V2i 3 
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V27 \V2 



cos ( -^=) sinh ( 

V2J \V2 



(C13) 



For a > 1/ y2 this leads to the following representation of 04(a) = #4,4 (a) (see Eq. (jClip ): 

poo 

g 4 (a) = 7T 3 / 2 / dt t 3 / 2 £ 4 ,4(-£ 4 )exp [-ta 2 ] 
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sin ^ 
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sin ( -arccotfv^a 2 — 1] j + cos ( -arccot["\/2a 2 — 1] 



1 + [V2a 2 - 1 

and the equation for the order parameter profile becomes 

f 

arctanh [<£*(*)] = -/3(A/i - AV(z)) + K I a 4 <j>*(z) + <" [<p*(z + 1) + <j>*{z - 1)] + A E 9^\ 9 ~ Z 'D^*( Z ') 



(C14) 



|z'-z|>2 



JC151 



where if = /?Jg r , and A = J'/J^.. In Eq. (|C 15|1 one has where 



a 4 = 7 + A( C4 -4), <" = 1 + A(cj n - -), (C16) 



C4 = 



F — 

4i 1 + n ■ 



4.900 



(C17) 
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and 



cf 1 = gi{l) + cYl~ 1.015 + 0.013 « 1.028. (C18) 



This shows that the contribution of cj n to df 1 is ca. 1%. 
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FIG. 14: The behavior of the function ffd(a) and its asymptote 
for d = a — 3 and d = a = 4. 

For the asymptotic behavior one finds 



/ \ " -i5 

g A (a —> 00) = —a 



33^ 
1024' 



,-13 



+ 0{a- 21 ). (C19) 



By setting g 4 = 0, Eq. (|C15[) reduces to the well-known 
form of the mean-field theory for systems with short- 
ranged fluid-fluid forces. (Actually in our case setting 
.94 = corresponds to a system with short-ranged in- 
teractions in z-direction and long-ranged ones within 
the planes perpendicular to z; only setting A = or, 
equivalently, J 1 = reduces the model to one with 
strictly short-ranged fluid-fluid interactions.) The stan- 
dard Ginzburg-Landau equation follows from taking the 
expansion arctanh(0 — > 0) = <f) + <f) a /3 + O(0 5 ). A contin- 
uum version of the equation follows from the replacement 
4>(z - 1) + 4>(z + 1) -> 2(j>(z) + (j)"{z). Obviously such a 
continuum version can be constructed also for a system 
with long-ranged forces by taking into account the func- 
tion 34(a) (which then has to be considered as a function 
of a continuous variable) . This exploits the fact that the 



function 34(a) is well defined for a > (provided that for 
one takes the principal values of arccot 
to be in the interval (0,7r)) and not only for a > 1, as 
actually needed for the lattice formulation of the theory. 
Therefore in the continuum formulation of the theory the 
integration can be extended to the full interval z £ [0, L\. 
This does not change the algebraic asymptotes of the 
density profiles. Thus in the continuum case the equa- 
tion for the density profile in the presence of long-ranged 
interactions is an integro-differential equation: 

0*(z) + i(0*(z)) 3 = i/3(A M -AV(z)) 



+K^a i( j)*(z) + a^ 



24,* (z) 



d 2 (j)*(z) 



dz 2 



+X / gi (\z~z'\)r(z')dz' 
Jo 



(C20) 



The behavior of the function g 4 (a) as well as its asymp- 
tote (Eq. <|C19P ) are shown in Fig.Q31 Equation (|C20j) 
has to be augmented by the boundary conditions 4>*(z — > 
0) = (f>*(z —> L) = 00. 

Note that, according to Eqs. <jC15|) (or (|C20|) in the 
continuum case) and (| A8[) , the only explicit L dependence 
in the behavior of the order parameter profile enters via 
the asymptotes of the functions V(z) and (74 (1^1) for \z\ ^> 
1. They will then lead to contributions of the order of 
L~ a in the behavior of p away from T c , i.e., there p ~ pb + 
0(L~ cr ). The last result implies that the entropy term 
in Eq. (|2.2p in this regime will produce contributions to 
Wex of the order of L~ a . Such contributions can then be 
neglected in comparison with the ones described by Ai 
(see Eq. CE]))> A $ ( see E Q- l|A10P ). and A^ s (see Eq. 
(|B3P ). Obviously this argument is not specific for d = 4 
but is generally valid. Finally, we repeat that in the case 
d = u = 3 also an analytical expression for gz(z) can be 
derived (Hi: 



53(a) 



V3 arctan 



V3 



In 1 - 
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(C21) 



The behavior of this function as well as its asymptote are 
also shown in Fig.[l4l 
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